Direct observation of Levy flight of holes in bulk n-InP 
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We study the photo-luminescence spectra excited at an edge side of n-InP slabs and observed from 
the broadside. In a moderately doped sample the intensity drops off as a power-law function of the 
distance from the excitation — up to several millimeters — with no change in the spectral shape. The 
hole distribution is described by a stationary Levy-flight process over more than two orders of 
magnitude in both the distance and hole concentration. For heavily-doped samples, the power 
law is truncated by free-carrier absorption. Our experiments are near-perfectly described by the 
Biberman-Holstein transport equation with parameters found from independent optical experiments. 
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The transport of minority carriers produced by optical 
excitation in semiconductors is usually well described by 
a diffusion equation, characterized by a diffusion length 
I 2 = Dt, with D being the diffusivity of carriers and r 
their lifetime. The diffusion process can be viewed as 
the result of a random walk in which every step Xi has 
the same probability density V(xi). In moderately doped 
direct-gap semiconductors, the transport of minority car- 
riers is mediated by the transport of interband photons 
with carriers recombining and re-emerging repeatedly in 
a process called the "photon recycling". The random- 
walk steps Xi are defined by the reabsorption length of 
photons produced in radiative recombination. If the sec- 
ond moment of the distribution (x 2 ) is finite, one has 
D oc (x 2 )/t,i, where r, is the mean time between steps 
(the radiative emission time r ra d)- This may lead to a 
substantial enhancement of D and I but the diffusion de- 
scription would still be applicable. 

If, however, the step distribution V(x) is heavy-tailed, 
asymptotically satisfying a power law 



P{x) ~ l/x 



7 + 1 



(1) 



with < 7 < 1, then the conventionally defined diffu- 
sivity diverges and the random walk is governed by rare 
but large steps pj. Such a transport, called the Levy 
flight, features an anomalously large spread in space and 
a "superdiffusive" temporal evolution [H-Q]. 

Superdiffusive transport of light has been known in 
the context of radiation trapping in propagation through 
media with narrow absorption lines. This phenomenon 
occurs in different systems, ranging from stars 0, 0] to 
dense atomic vapors 0, such as gas lasers, discharges 
and hot plasmas. It has been studied for many decades, 
starting from the theoretical papers by Biberman and 
Holstein Recently, Levy flights of photons were di- 
rectly observed in engineered optical materials floj | and 
in hot vapors jllj . Direct-gap semiconductors represent 
a new and exciting "lab" system for studying Levy flight 
of photons and/or minority carriers 

mill- 

We have investigated [12[ the step distribution V{x) 
for the photon-assisted transport of holes in moderately 
doped n-InP and found that it asymptotically obeys the 



power law (QJ with 7 = 0.7±0.I. For heavier doping, the 
power law is truncated by free-carrier absorption, 



V{x) ~ (l/x 7+1 ) e~ a f 



(2) 



where a/ c [cm ps 1.3 x 10 
absorption coefficient in InP 



~ 18 iV,v is the free carrier 



mm 



The resultant random walk is limited by the loss of 
carriers in nonradiative recombination (of rate t^. 1 ) . The 
relative rates of recombination are characterized by the 
radiative efficiency r\ = T nr /(r nr + r ra( j) and the typical 
number of steps by the recycling factor <I> = 77/(1 — rf). 
Due its high radiative efficiency (with $ reaching 10 2 for 
moderately doped samples), the n-doped InP is ideally 
suited for studying the Levy flight of holes 12-[l5j]. The 
key parameters (7 and can be controlled by vary- 
ing the doping Nd and the temperature. The truncation 
length aj^ can also be controlled by varying Nd- 

In earlier photoluminescence experiments with n-InP 
we found an indirect evidence of Levy flight of holes by 
analyzing the ratio of transmitted and reflected lumines- 
cence spectra across a thin flat wafer 



12] 



A very different geometry is used in this work: the lu- 
minescence is excited at an edge side of the wafer and 
observed from the broadside, Fig. [TJ The luminescence 
intensity measured as a function of the distance x from 
the edge is directly proportional to the hole distribution 
p{x) in the sample. In a moderately doped sample, we ob- 
serve a power-law distribution, providing an unambigu- 
ous manifestation of the Levy flight of holes. In sam- 
ples of heavier doping, p{x) evolves toward a distribution 
characteristic of truncated Levy flight. 

We have studied three n-InP samples [l6| of thickness 
d=350 fxm and doping N d = 0.3, 2 and 6xl0 18 cm -3 
(samples I, II and III, respectively). A 808-nm laser 
beam was focused (see Fig. [T]) on the 7-mm edge side 
of the sample using a cylindrical lens. The excitation 
flux was uniform along the edge. Excitation photon en- 
ergy E = 1.53 eV (above the interband absorption edge 
at Eg w 1.35 eV) ensured the initial hole generation in a 
thin layer near the surface. Luminescence emitted from 
the broadside was captured by a lens to project the image 
on a CCD camera. A razor blade was installed near the 
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FIG. 1: Geometry of our photoluminescence experiment. To 
avoid scattering by the edges, the laser beam is focused on a 
narrow strip zo ± Sz in the 7 mm side of the sample, where 
Sz ~ 50/im and zq ~ lOO^m, counting from the top surface. 



excitation edge to obscure the luminescence emitted in 
the direction of the camera from the edge surface itself. 

We have also measured the luminescence spectra at dif- 
ferent distances x from the excitation. These normalized 
spectra are shown in Fig. [5]for sample F In the range of x 
from 0.5 to 3 mm, the spectra remain essentially the same 
and in agreement with the calculated spectrum (dashed 
line) for a near- uniform hole distribution p(z) = p(xq, z) 
» const, as expected for any fixed xq > d (see below). 

The absorption spectrum a{E) of sample I is shown 
in Fig. [5] by a dash-dotted line. It exhibits an exponen- 
tial Urbach behavior, extending down to a rj 0.5 cm . 
At lower E, the Urbach tail, characteristic of interband 
absorption m, is masked by the free-carrier absorption, 
a(E) = a.i{E) + ctf c . At high E, the Urbach exponent 
saturates above E g . With the higher doping, the absorp- 
tion edge is shifted to higher energies due to the Moss- 
Burstein effect. It is notable that the observed lumi- 
nescence spectrum lies fully in the Urbach tail region of 
the absorption spectrum — in contrast with the intrin- 
sic emission spectrum described by the quasi-equilibrium 
van Roosbroek-Shockley [TtJ relation (VRS), 

S VRS (E) ~ a t E s e~ E ' kT , (3) 

which has a maximum above E g . 

The luminescence intensity distributions I(x) were ob- 
tained by scanning the CCD image near its center along 
a line parallel to x. To reduce random fluctuations, the 
distributions were averaged over 20 scans. Additionally, 
the CCD camera was shifted both in the y and x direc- 
tions, with the corresponding scans averaged again to re- 
duce irregularities in the response of different pixels. The 
observed I(x) was strictly proportional to the excitation 
laser power and hence to the hole concentration p{x). To 
stay within the linear-response range of CCD over the 
entire range of x, we first used a neutral filter to reduce 



I(x) near the excitation edge and then the razor blade 
was shifted along x to obscure the brightest parts of the 
slab just near the edge. Subsequently, the filter was re- 
moved to get a measurable signal far from the edge. The 
residual dark noise of the camera was subtracted. The 
resulting distributions for all three samples are presented 
in Fig. [3] 

For sample I at distances x > 0.5 mm, a power law 
I(x) ~ l/x 1+1 is clearly observed. This is most easily 
seen on the log-log scale in the inset of Fig. [3] The best 
fit for 7 = 0.7 ± 0.1 agrees with the index 7 obtained 
earlier [12[ for V(x). The power law is in clear contrast 
to an exponential decay I(x) ~ exp(— x/l) expected for a 
(photon-assisted) diffusive spread of holes, even account- 
ing for any enhancement of the diffusion length I by recy- 
cling (a diffusive curve for /=210 /im is shown in Fig. [3] 
by the dashed line). 

For samples II and III with higher doping, the heavy 
tails are also clearly seen. However, the power-law dis- 
tribution is truncated at large distances. This effect cor- 
relates with the increasing af c (Nd). 

The nature of the emission spectra is discussed in 
detail in the review [l3j]. The room-temperature in- 
trinsic emission spectrum in n-InP is well described by 
Eq. ([3]). The observed spectrum S(E) is modified (fil- 
tered) by rcabsorption on its way out of the sample, 
S(E) = S V rs x F(E), where F(E) is the radiation fil- 
tering function F(E) = F\(E) T{E) which depends on 
the hole distribution p{z) across the wafer and is affected 
by reflections from the sample surfaces. The one-pass 
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FIG. 2: Luminescence spectra observed from the top side 
of sample I at varying distances from the edge (solid lines). 
Dashed line shows the calculated spectra assuming an intrin- 
sic van Roosbroek-Shockley emission spectrum, modified by 
reabsorption of radiation in the sample. Also shown is the ab- 
sorption spectrum a(E) exhibiting an Urbach tail of Ui below 
E g in a wide range of variation. The Urbach-tail dependence 
is masked by the free-carrier absorption only for E < 1.26 eV. 
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FIG. 3: Distribution of the luminescence intensity I(x) and 
hole concentration p(x) oc 7(x) for three differently doped 
samples (dots). Results of calculations using optical parame- 
ters of the samples are shown by solid lines. For comparison 
also shown (by a dashed line) an exponential distribution with 
enlarged diffusion length I — 210 fjja. For the low-doped sam- 
ple I, the distribution is a power law, with the index 7 readily 
obtained from the slope of the log-log scaled graph in the in- 
set. For samples II and III with higher Nd, the distribution 
is truncated by free-carrier absorption. 



filtering function F\(E) is given by mm 



Fx{E)= / p(z)exp[-a(E)z]dz 



(4) 



The factor T(E) = [1 - Rex^{-a{E)d)]- 1 accounts for 
the multiple surface reflections. Note that due to the high 
index contrast the outgoing radiation propagates close to 
the normal direction to the surface and the surface reflec- 
tion coefficient R(E) w 0.33 varies little in the narrow 
energy range of the emission line. Multiple reflections 
are noticeable in the red wing of the spectrum. 

One can expect that for x > d the distribution p(z) 
is uniform across the sample, except for small regions 
near the surfaces, where it is distorted by surface re- 
combination. This assumption is confirmed by solving 
the one-dimensional diffusion equation with a recycling 
term [see Eq. ([5]) below] for p(z) with uniform genera- 
tion at a given x. The calculated emission spectrum is 
fully determined by parameters of the absorption spec- 
trum. Results of the calculations, shown in Fig. [5] by a 
dashed line, agree with experiment near perfectly, with- 
out adjustable parameters. Position of the emission line 
maximum can be calculated analytically by substituting 
p(z) = const into Eq. (jl]). Near the center of the emis- 
sion line the absorption coefficient is described by a sim- 
ple exponent a = aoexp[(E — E g )/A] with the param- 



mcV, known from optical studies [FJ, |18( . The emission 
line peak position E mSLX is then found from the equation 
dS{E)/dE = 0: 



A\n(a d/s) 



(5) 



where s is a non-zero solution of transcendental equation 
(kT/A)s = [exp(s) - 1]. This gives £ max = 1.303 eV, 
in agreement with experiment. This agreement further 
supports our assumption of the shape (j3]) for the intrinsic 
spectrum and confirms the proportionality p(x) oc I(x). 

Next we discuss the hole distribution p(x) along the 
sample. Since the excitation is restricted to a narrow re- 
gion near the edge surface, and all surfaces are highly re- 
flective (with a narrow radiation escape cone), this distri- 
bution can be found from the one-dimensional Biberman- 
Holstein stationary transport equation [H], HI, 0, [l3| 



- D 



d 2 p(x) p{x) 
dx 2 T 



G{x) + 1 

T 



p{x')V{\x-x'\)dz' 
(6) 

where D is the ordinary hole diffusivity and r 1 is the 
total recombination rate. The last term in the right-hand 
side accounts for the photon recycling. The reflection 
from the edge face at x = is included by assuming a 
symmetric distribution p(x) = p(—x) and extending the 
integration to —00 (which supplies an image source for 
every radiative recombination event in the sample). Here 
V(\x~x'\) is the probability, averaged over the plane x = 
const, of a hole to generate another hole at a distance 
\x — x'\ by the emission- reabsorption process (see [IH): 



V(\x\) 



Af(E)ai{E)Ei(l,a(E)\x\)dE , (7) 



where Ei(l,x) = J dt V 1 exp(— xt). The integrand in 
Eq. ((?]) is a product of probabilities for (i) emission of a 
photon at energy E, described by a normalized spectral 
function J\f(E) oc E~ x Svrs(E), (ii) propagation of this 
photon from a point at distance x to another point at 
distance x', described by the factor Ei(l, a(£ , )|x|), and 
(iii) interband absorption of this photon, represented by 
the factor 0*4 (E) . Hence, neglecting the very small effects 
of hole displacement in their thermal motion, the func- 
tion V(\x — x'\) is the single-step length distribution for 
holes. The probability Pdxl) given by Eq. (J7J) is again 
fully determined by the absorption spectrum and can be 
calculated numerically 12|,ll3|. For all x > 1/im, it obeys 
the power law ([2]) with an exponential factor allowing for 
free-carrier absorption. The index of the distribution is 
fully determined by the Urbach-tail part of the spectrum 
and is given by 7 = 1 — A/ kT. 

Solution of Eq. ([6]) can be obtained by a Fourier trans- 
formation. For G = Gq5(x) it is of the form 



p(x) 



G t 



cos(fcx) 



1 + l 2 k 2 ~ r]F(k) 



dk 



(8) 
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Here I = \J~Dt is the ordinary hole diffusion length, and 
F(k) is the Fourier-transform of P(|x|). For x ^s> I the 
distribution p(x) is fully determined by F(k) (and hence 
in turn by the absorption spectrum). The only additional 
parameter rj has been already determined for each sample 
- from the time-resolved luminescence kinetics experi- 
ments [lit and, independently, by analyzing the ratio of 
transmission-to-reflection luminescence spectra 12 1 . 

An alternative to solving Eq. © is to do Monte Carlo 
modeling of the hole distribution using the single-step 
probability j7j [L2J, [l3| . The result coincides with flSJ at 
distances x < 500 yum. However, at considerably larger 
distances of interest to us here, the Monte Carlo approach 
becomes noisy and therefore less reliable flit ]. 

Results of the numerical calculations of p(z) using ([8]) 
for all samples are shown by the solid lines in Fig. [3] 
demonstrating an excellent agreement with the experi- 
mental data, except for a slight discrepancy for the heav- 
iest doped sample III (see below). 

Main features of the distribution are revealed by an 
analytic approximation that can be derived directly from 
Eq. ([6|). At large distances, due to the heavy tail of 
■p(|a;|), one can solve Eq. (JB|) by sequential iterations — 
making in all terms of the resultant series the "longest 
step approximation" [f| (i.e. choosing one of the steps 
equal the total distance x). This gives 



p{x) 



<f>V(x) 



[1 + $V e (x)]< 



(9) 



where V e [x) = J V(x')dx' is the probability of escape 
beyond x in one step. Equation (|9|) provides a good 
approximation to the exact solution in the entire range 
of x. It works with or without truncation of the Levy 
flight by free-carrier absorption. Similarly to the "sta- 
ble distribution" [l|, distribution © asymptotically re- 
produces the one-step probability (enhanced by the recy- 
cling factor <&) and is modified at smaller distances. Since 
$P e (x) is the escape probability in $ attempts, condi- 
tion $>T> e (xf) = 1 gives the distance Xf to the "front" 
of p(x) — beyond which the holes appear predominantly 
in one step. The front Xf is clearly seen in the inset of 
Fig. |3l as the point of maximum curvature on the log- 
log plots for all samples. The observed Xf « 200 /jm 
manifests the large values of $ in these samples. For 
x ■< Xf, Eq. © correctly predicts the short-distance 
asymptote [13[, p(x) ~ 1/.t 1-7 but the corresponding 
data are obscured in our experiment by deviations from 
the one-dimensional geometry (20j . For sample III, the 
effective front distance Xf is somewhat larger than could 
be expected from earlier experiments. 

In conclusion, we have studied the stationary hole dis- 
tribution p{x) over distances x of several mm from initial 
photoexcitation. We observe a heavy-tailed decline of 
the luminescence intensity with no change in the spec- 
tral shape. For a low-doped sample, we find a power-law 



distribution p(x) precisely accounted for by the Levy- 
flight transport of holes mediated by photon recycling. 
In heavier-doped samples, the power law is truncated by 
free-carrier absorption, which manifest itself only at dis- 
tances x ~ a~jl corresponding to the free-carrier absorp- 
tion. The anomalous transport should be important in 
all semiconductor crystals with high radiative efficiency 
and may have practical implications in optoelectronic de- 
vices [13| . It is of interest to study enhancement of this 
effect at lower temperatures. 
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